All-electron self-consistent GW approximation: Application to Si, MnO, and NiO 
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We present a new kind self-consistent GW approximation (scGVF) based on the all-electron, full- 
potential LMTO method. By iterating the eigenfunctions of the GW Hamiltonian, self-consistency 
in both the charge density and the quasiparticle spectrum is achieved. We explain why this form 
of self-consistency should be preferred to the conventional one. Then some results for Si are shown 
as a representative semiconductor, to establish agreement with a prior scGW calculation. Finally 
we consider many details in the electronic structure of the antiferromagnetic insulators MnO and 
NiO. Excellent agreement with experiment is shown for many properties, suggesting that a Landau 
quasiparticle (energy band) picture of MnO and NiO provides a reasonable description of electronic 
structure even in these correlated materials. 

PACS numbers: 71.15.-m,71.15.Ap,71.20.-b,75.50.Ee 



The GW approximation (GWA) of Hedin[l| is gener- 
ally believed to accurately predict excited-state proper- 
ties, and in particular improve on the local density ap- 
proximation (LDA), whose limitations are well known, 
e.g. to underestimate bandgaps semiconductors and in- 
sulators. Usually GWA is computed as 1-shot calculation 
starting from the LDA eigenfunctions and eigenvalues; 
the self-energy S is approximated as S = iG^^^W^^^ , 
where G^^^ is a bare Green function constructed from 
LDA eigenfunctions, and ly'"'^^ is the screened Coulomb 
interaction constructed from G^^^ in the random phase 
approximation (RPA). However, establishing the validity 
of the 1-shot approach has been seriously hampered by 
the fact that nearly all calculations to date make further 
approximations, e.g. computing S from valence electrons 
only; the plasmon-pole approximations; and the pseu- 
dopotential (PP) approximation to deal with the core. 
Only recently when reliable all-electron implementations 
have begun to appear, has it been shown that the 1-shot 
GWA with PP leads to systematic errors [3, 0, There 
is general agreement among the all-electron calculations 
(see Table IJ) that the F-X transition in Si is underes- 
timated when E = iG^^^W^^^ . And we have shown 
prcviously0 that the tendency for E = iG^DA^LDA 
underestimate gaps is almost universal in semiconduc- 
tors. This is reasonable because W^^^ overestimates 
the screening owing to the LDA small band gaps. G con- 
structed from quasiparticles (QP) with a wider gap (e.g. 
a self-consistent G) reduces the screening, and therefore 
generates GW with a wider gap. 

However, there are many possible ways to achieve self- 
consistency. The theoretically simplest (and internally 
consistent) is the fully self-consistent scheme [scGW), 
which is derived through the Luttinger-Ward functional 
with the exchange-correlation energy approximated as 
the sum of RPA ring diagrams. Then W is evaluated 
as W ^ v{l — vP)~^ with the proper part of the polar- 



ization fmiction P= —iG x G. However, such a construc- 
tion may not give reasonable W , resulting in a poor G, 
for the following reason. If S is w-dependent, G can be 
partitioned into a QP part and a residual satellite part. 
The QP part consists of terms whose energy-dependence 
varies as Zi/{u!— eiiiVi), where e^, F^, and Zi are respec- 
tively the QP energies, inverse lifetimes, and renormaliza- 
tion factors {Zi<l; typically between 0.7 and 1). The QP 
parts are thus weighted by factors Z; the residual weights 
1 — ^ go into the plasmon-related satellite parts, high in 
energy. Thus P — —iGxG contains contributions from the 
particle-hole pair excitations as does —iG^^^xG^^^ , but 
reduced by the products of two Z factors, one from oc- 
cupied and the other from unoccupied states. However, 
this construction of P is not consistent with Landau's 
quasi-particle theory, which insists that one-particle ex- 
citations remain meaningful (at least near the Fermi en- 
ergy). Based on the theory, we should instead evaluate 
the QP contributions to P without Z factors, as they 
dominate the static screening W{uj^O). Inclusion of Z 
can lead to W{0) being underscreened; moreover W{lu) 
does not satisfy the / sum rule Consequently W{0) 
will be overestimated resulting in a tendency to overesti- 
mate band widths, as is well known in the extreme case 
(Hartree-Fock), where there is no screening of W. 

Indeed Holm found that sc G IF overestimates the band 
width in the homogeneous electron gas|^. Very recently 
Ku presented a sc GIF calculation which similarly overes- 
timates the valence band width in Si and Ge; see Table HI 
0. Another practical justification for the argument that 
a bare G without Z should be used when we construct P 
as —iG X G, is that W^^^ is already known to be rather 
good if we add some enlargement of band gap by hand 
to correct for errors in the LDA 7J. 

For this reason, we do not adopt the full sc GIF scheme, 
but construct two kinds of constrained self-consistent 
GIF methods. For a set of trial eigenfunctions and 
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quasi-particle energies {i/^qm^qn}, we can calculate the 
one-particle Green function G, from which we can cal- 
culate self-energy E^^, (w) in the expansion of {ipqn} in 
the GWA. Then we generate an energy-independent and 
hermitian S^,^, in one of two ways: 



Re 



(Cqn) + (eq„') /2 



(2) 



where Re means that we take only the hermitian parts. 
With this S^„/, we can construct a new density n(r) 
and corresponding Hartree potential, and proceed to a 
new set of {tpqn, eqn}- The procedure can be iterated to 
self-consistency. Our method is not related to the LDA 
(though in practice the LDA is used to make a starting 
guess for E and the augmented- wave basis set). These 
two schemes differ in the treatment of the off-diagonal 
parts of S„„/ ; but both restrict the potential to be non- 
local, hermitian and w-independent. Thus the problem in 
the full scGFFis avoided; also the numerical computation 
becomes rather stable. In Ref.^j], off-diagonal matrix ele- 
ments of E„„, were completely neglected. However, these 
are important for MnO and NiO: eigenfunctions and the 
density can not be changed from LDA if we neglect them. 
We find that converged QP energies differ in these two 
schemes by small amounts (less then < 0.02 eV for Si, 
typically ~ 0.1 eV for NiO), which is within the resolu- 
tion of the method (~0.1 eV). 

Our implementation is based on the method of Ref. (2j . 
W is expanded in a mixed basis which consists of two 
contributions, local atom-centered functions (product ba- 
sis) confined to muffin-tin spheres, and plane waves with 
the overlap to the local functions projected out. The 
former can include any of the core states: thus the va- 
lence and core states can be treated on an equal foot- 
ing and the contribution of the latter to S included. 
We calculate the full energy dependence of W without 
the plasmon-pole approximation. This approach shares 
some features in common with both the full-potential, all- 
electron plane- wave based methods 0,0 and the product- 
basis method^, combining the advantages of each, e.g. 
efficient treatment of localized valence electrons. 

Results for Si are shown in Table |l| Agreement be- 
tween the three all-electron methods is generally excel- 
lent. The G^^^W^^^ gaps are ~0.3 eV smaller than 
experiment; the scGFFgaps fall much closer. As we will 
show elsewhere, most properties of weakly correlated sys- 
tems calculated with the present scGW method (funda- 
mental and higher-lying gaps, valence bandwidths, ef- 
fective mass, position of deep d levels) are in excellent 
agreement with experiment, with small but systematic 
residual errors. 

Turning to the TM oxides, we first consider MnO be- 
cause it is less correlated. Fig.ncompares the scGWen- 
ergy bands and corresponding DOS to the LDA and the 
(^LDA|^LDA gj^p^ rpj^g conduction band at F is evidently 



TABLE L Minimum energy gap Eg and selected energy eigen- 
values for Si, relative to T'2Sv (eV). Three all-electron methods 
are shown: linearized augmented-plane-wave (LAPW) and 
projector-augmented-wave (PAW) approaches, and LMTO 
(this work). The PAW calculation included valence electrons 
only. The last row compares the Ge valence bandwidth. The 
results of this work differ slightly from Ref. Q because a large 
basis set (50 orbitals/atom) was employed in the present work. 



PAW[3] 
(Giy)LDA 


LAPW|^ This work 
(GVF)^°^ scGW {GW)^°^ scGW 


Exp. 


Eg 0.92 


0.85 


1.03 


0.92 


1.14 


1.17 


Xi, 1.01 






1.06 


1.30 


1.32 


Lie 2.05 






2.00 


2.26 


2.04 


Fisc 3.09 


3.12 


3.48 


3.11 


3.40 


3.38 


Flu 


-12.1 


-13.5 


-12.1 


-12.3 


-12.5 


Fi4Ge) 


-13.1 


-14.8 


-12.9 


-13.1 


-12.6 



a dispersive band of sp character. Above this, fall the t2g 
bands (~6-9 eV); still higher at ~10 eV is a narrow eg 
band, whose width is ~3 eV. Thus, the itinerant and d 
bands are well separated. The minimum gap is 3.5 eV, in 
good agreement with the BIS gap0 (3.9±0.4 eV). The 
BIS spectrum also shows a peak at ~6.8 eV, which prob- 
ably corresponds to a convolution of the peaks of t2g 
symmetry seen in the DOS at 6.6 eV and 7.3 eV. These 
bands are in stark contrast to the LDA, which shows the 
t2g and Cg bands overlapping and hybridizing with the 
sp band at 1 to 4 eV. 

LDA and scGIF valence bands are more similar 
(Fig. . In the LDA there is a narrow upper Cg band 
at 0.1 eV below the valence band maximum (VBM), and 
another one at VBM— 5 eV. Both weakly hybridize with 
the O 2p band. The scGIF pushes the upper eg band 
down to VBM- 0.5 cV, so that the VBM takes more O 
2p character, and the band at VBM— 5 eV takes more 
Mn d character. The splitting A^, between the upper 
eg level and the t2g level widens from 1.0 eV(LDA) to 
1.7 eV(scGlF), in good agreement a photoemission mea- 
surement of 1.9 eV|9|. An approximately similar picture 
emerges from a model GW calculation of Massidda|lo|. 
the most important difference being that the model GW 
d conduction bands fall ~1 eV lower than ours. 

Fig. 121 compares the sc GIF energy bands for NiO along 
the [110] and [100] lines to ARPES data of Shen (ll| for 
the valence bands, and to the LDA conduction bands. 
The right panel shows the density-of-states (DOS) for 
both LDA and scGIF , and Fig.Olshows the total DOS 
resolved into components. Also shown in the top panel 
are BIS data taken from Ref. fl^ 

Several features are of interest: 

1. The conduction-band minimum falls at the F point; 
the VBM falls at the point (1/2,1/2,1/2) (not shown). 
The calculated minimum gap is 4.8 eV. 

2. The sc GIF conduction bands are a mixture of a disper- 
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FIG. 1: scGW energy bands and DOS of MnO. Solid lines: 
scGW bands; dashed lines: LDA bands. The VBM is set to 
energy 0. Circle and triangle at T: BIS and G^°^W^°^ gaps. 
Right panel shows the corresponding DOS. Peaks at -0.5 eV 
(-0.1 eV in LDA) and -5 eV are the nearly dispersionless Cg 
bands. Peaks at -2.2 eV (-1.2 eV) and 6.6 and 7.3 eV (1.7 
and 1.9 eV) derive from Mn t2g states. 





(1,0,0) 



FIG. 2: scGW energy bands and DOS of NiO. Valence-band 
maximum is set at energy 0. Solid lines: scGVF bands; dashed 
lines: LDA bands (only conduction bands are shown). Circles 
show photoemission data of Ref. lid . Right panel shows the 
corresponding total DOS. 



sive band composed of sp approximately equally weighted 
on the Ni and O sites, and a nearly dispersionless Cg state 
(see discussion of EELS, below). Peaks in the BIS spec- 
trum labeled '1', '2' and '3' closely coincide to those in the 
scGFF total DOS, apart from a constant shift of 0.8eV. 

3. The sc G valence bands are in very good agreement 
with experiment: indeed they agree as well with the Shen 
data as the latter agrees with an independent experiment 
by Kuhlenbeck et. al.^^ (not shown). 

4. There is an increased dispersion in the valence bands 
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FIG. 3: DOS and EELS of NiO. Solid lines: scGW data; 
dashed lines: LDA data. Top panel: total DOS, together 
with BIS data of Ref. [12 (circles). Panels 2 and 3 show the 
Ni t2g and Cg partial DOS, with positive DOS showing major- 
ity spin and negative showing minority spin. Panel 4 shows 
the O sp partial DOS; panel 5 compares the calculated and 
measured |l^ EELS spectra from the O Is level. 



relative to the LDA at the VBM because the nearly dis- 
persionless Ni t2g levels are pushed down. Thus the VBM 
acquires somewhat more O 2p character. This supports 
the generally accepted picture that the LDA too heavily 
favors the Mott-Hubbard picture. 

The bottom panel of Fig. O compares the EELS spec- 
trum calculated as described in Ref. Iisl computing the 
excitation from the O Is core. Calculated data were con- 
volved by a Gaussian of 0.5 eV width (which was the 
resolution reported) to compare to experimental data by 
reported by Dudarev et. al.|jl^. The calculated results 
were shifted to align the spectra with the DOS. For pur- 
poses of comparison, Dudarev's data was shifted by 526 
eV to align the peaks with the scGiy results. (Peaks la- 
beled '1', '2' and '3' should correspond to the peaks with 
the same labels in the BIS spectrum; indeed the EELS 
peaks and BIS peaks almost perfectly align if the EELS 
data is further shifted by 0.8 eV.) 

Apart from the 0.8 eV shift, the EELS data is in ex- 
cellent agreement with the sc G results. Spacings be- 
tween the three peaks agree to within ^0.1 eV, and the 
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spectral weight under each peak (estimated by numerical 
integration) also agree closely. This establishes that the 
scGFF relative positions of the sp and Ni d bands are cor- 
rectly predicted. This is a significant result, because the 
relative positions of the sp and d bands is a rather a del- 
icate quantity^]. In contrast, the LDA underestimates 
the spacings between peaks 1 and 2 by '^1.5 eV and be- 
tween peaks 2 and 3, by ~0.7. Moreover, it overestimates 
the spectral weight of the first peak by a factor of ^2. 
This result is also significant, because the EELS spectra 
largely reflect the O 2p partial DOS. Without coupHng 
between the Ni eg level, the itinerant band would adopt 
a simple parabolic form; thus the amplitude of first peak 
is a reflection of the hybridization between the Ni Cg and 
the itinerant band. The fact that scGW^gets the correct 
weight for this peak establishes that it accurately esti- 
mates this coupling, while the LDA overestimates it by 
a factor of 2. 

Many of the results found here confirm many conclu- 
sions drawn in a model calculation (l^l, as well as 
various LDA-Kf7 calculations [H El Ei, 

both of which 

may be viewed as model approaches to the present the- 
ory. Some significant differences do arise. The rela- 
tive positions of different bands and the energy gaps de- 
pend rather sensitively on the choice of parameters in 
the model approaches. For example, in Ref. E3i 
0-derived sp conduction-band appears to fall at 2.8 eV, 
~2 eV below the middle of the Cg level when U is assumed 
to be 5 (somewhat lower than the constrained LDA es- 
timate, eV). Massidda's model GW calculation^^ 
shows the sp band ~1 eV above the eg. 

To what extent does the Landau QP picture based on 
the preceding scGW^ results fail to describe the true elec- 
tronic structure of MnO and NiO? We have shown that 
a great deal is correctly described, including many de- 
tails of the valence and conduction bands. The main 
discrepancy is with XPS measurements. For optics, 
the peak in Im(e) corresponding to the gap in NiO is 
about VBM-l-4 eV whereas this peak is at VBM-l-5 eV 
in the BIS data. But Im(e) is directly related to the ex- 
citonic process or correlational motion of electron-hole 
pair. Such a correlation can shift Im(e) downward. The 
difference between the two experiments can be due to 
this correlation. So the poles of the true Green's func- 
tion (which are reflected in the DOS) should correspond 
to the unoccupied d position at VBM+5 eV as is shown 
in BIS. Peak '1' in the scGFFDOS falls slightly higher 
than experiment, at ^5.8 eV. If we include correlation 
beyond RPA, e.g. inclusion of ladder diagrams, the band 
gap in the Green's function may be reduced about 1 eV, 
as estimated by Takahashi and Igarashi20j. Thus, it 
would seem that the RPA explains quite well the impor- 
tant experimental data, apart from a slight tendency to 
underestimate screening of I^jlJl. We have not yet at- 
tempted to include excitonic effects, so we cannot say to 
what extent photoemission data can be explained within 



TABLE II: Magnetic moments and minimum gaps in MnO 
and NiO. 





moment 


bandgap 


Compound 


LDA scGW^ Expt 


LDA scGW Expt 


MnO 


4.48 4.76 4.6 


0.78 3.5 3.9±0.4 


NiO 


1.28 1.72 1.9 


0.45 4.8 4.3 



the RPA, though estimates in a model context were rea- 
sonably successfulE^. Thus, we believe that the band 
pictureji^l for NiO is a reasonable starting point for the 
description of the electronic structure of NiO, much bet- 
ter than previously thought, and in many respects more 
appropriate than the ligand-field picture. 

This work was supported by the Office of Naval Re- 
search. 
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